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The complete spectrum of plasmons of the two interpenetrating plasma streams is found in a 
closed analytic form. The orientation of the wave vector with respect to the stream direction is 
arbitrary and the plasmas, which are assumed to be collisionless and spatially homogeneous, can be 
nonrelativistic, relativistic or even ultrarelativistic. Our results apply to the electromagnetic plasma 
of electrons and passive ions and to the quark-gluon plasma governed by QCD. 


I. INTRODUCTION 


A plasma system, which is a complex of particles and fields, supports a wide variety of waves. Using quantum 
terminology one says that the plasmas reveal rich spectra of collective excitations or modes. The spectra carry 
information about the thermodynamic and transport properties of equilibrium plasmas and about the temporal 
evolution of non-equilibrium systems. The waves or modes corresponding to oscillations of charge density or currents 
are usually of the highest frequency and thus they play a particularly important role in the plasma dynamics. In the 
case of electromagnetic plasmas we deal with electric charges and currents and the electromagnetic waves are classical 
representation of quasiphotons. In the quark-gluon plasma governed by QCD there are color charges and color currents 
and the quasigluon is the analog of the quasiphoton. We call the collective modes related to the oscillations of charges 
and currents as plasmons, not limiting this term to the equilibrium plasma. 

Plasmons, which are found as solutions of dispersion equations, have been extensively studied for various plasma 
configurations over decades and there is a huge literature about the problem, see e.g. the textbooks mia. However, 
there is a rather limited number of exact analytic solutions of dispersion equations and there are only exceptional 
cases when the complete spectrum of plasmons can be found in a closed analytic form. We present here a general 
and exact solution of the dispersion equation of the two-stream system for an arbitrary orientation of the wave 
vector with respect to the stream direction. The plasma system is collisionless and spatially homogeneous, but it can 
be nonrelativistic, relativistic or even ultrarelativistic. Our results apply to both electromagnetic and quark-gluon 
plasmas. The complete spectrum, which is given in a closed analytic form, consists of eight modes (four pairs of modes 
of opposite sign) which are either real or imaginary. For any orientation of the wave vector, there is one unstable 
mode. When the wave vector is perpendicular to the stream direction, we have the magnetic filamentation mode 
which continuously changes into the electrostatic longitudinal mode when the wave vector becomes parallel to the 
stream. 


Two real modes exhibit the phenomenon of mode coupling, see e.g. the handbook [3], when the wave vector is 
almost perpendicular to the stream or the stream velocity approaches the speed of light. Then, the two dispersion 
curves approach each other but they do not cross. The effect is encoded in the solutions of the dispersion equation we 
found and to reveal it no additional reasoning or approximate methods, which are typically used [3], are not needed. 
The whole spectrum is rather rich and nontrivial and it changes qualitatively with the orientation of the wave vector. 

The two-stream plasma system, which occurs in various plasma physics experiments, is repeatedly discussed in the 
literature, see e.g. the textbooks mm- In context of the quark-gluon plasma it was first studied in mm- However, 
as far as we know, the complete spectrum of plasmons in such a system has not been found in closed analytic form 
and thus our results are new and original. And since analytic form of the spectrum greatly simplifies derivation of 
various plasma characteristics, the results are also practically useful. 


Presentation of our analysis is organized as follows. In the forthcoming Sec. |TT] we first formulate the problem by 
writing down the general dispersion equation of plasmons. Then, the momentum distribution of plasma constituents, 
which controls the dielectric tensor, is briefly discussed and an explicit expression of the tensor, which enters the 
dispersion equation, is derived. In Sec. m we introduce the method to solve the dispersion equation which is of the 
form det[£] = 0. Instead of computing the determinant we rather invert the matrix £ and look for poles of £ _1 . 
An analysis of the collective modes starts in Sec. IV with the discussion of two special cases. Then, we derive the 
complete set of exact solutions of general dispersion equations in Sec. [V] The spectrum of plasmons is found for any 
orientation of the wave vector. Sec. |VI| discusses the case, which is interesting both physically and mathematically, 
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when the stream velocity approaches the speed of light, 
remarks. 


In Sec. VII we summarize our study and make some final 


Our analysis of the two-stream system is methodologically very close to the extensive study of plasmons in the 
system where the momentum distribution of plasma constituents is obtained from the isotropic one by stretching or 
squeezing it in one direction [B], Since the two-stream distribution does not belong to this category, the technique 
to invert the matrix E needs to be modified. However, there are some repetitions with the article [B] which are 
unavoidable to make the present paper self-contained. We use notation very similar to that in Ref. [6] with the 
natural units, where H = c = 1, and the indices i,j , k = 1, 2,3 which label the Cartesian spatial coordinates. Lengths 
of vectors like k or u are denoted as k and u. 


II. FORMULATION OF THE PROBLEM 


The linearized chromodynamic or Maxwell equations of the Fourier transformed (chromo-)electric field can be 
written in the form 


E lJ (w,k)F; J (w,k) =0, 


(1) 


where the matrix E is defined as 


E y (w, k) = -k 2 6 ij + EkP + wV-V, k), 


( 2 ) 


ui is the frequency, k denotes the wave vector and e lJ (ui, k) is the (chromo-)dielectric tensor, see e.g. the review article 
[7J. External charges and currents are absent in the system under study and color indices, if needed, are suppressed 
in Eqs. ([TJ [2]) . A solution of Eq. (jT]) exists if 

det[E(w, k)] = 0, (3) 

which is the general dispersion equation. Its solutions w(k), which represent plasmons, are, in general, complex but 
the wave vector k is assumed to be real. There are the transverse plasmons, for which the electric field is transverse 
to the wave vector k, and longitudinal plasmons with the electric field parallel to k. The transverse modes correspond 
to oscillations of current, and the longitudinal ones to oscillations of charge density. A mode is called unstable if 
3w(k) > 0, because the amplitude ~ e 3W k )t grows exponentially in time. When 3w(k) < 0, the mode is stable 
and it is damped if 3w(k) < 0, as its amplitude decays exponentially in time. The mode is called overdamped , when 
additionally it is pure imaginary. 

To solve the dispersion equation, the dielectric tensor is needed. As discussed in e.g. [7 , the tensor for a locally 
chargeless anisotropic plasma in the collisionless limit equals 


e y (w,k) = S ij + 


2w 


d 3 p v l 

(27r) 3 ui — v • k -(- i0+ 



kw, • , 
- )<P fc + 


v^k k 

to 


v£/( p), 


( 4 ) 


where p, E p and v = p /E p are the momentum, energy and velocity of plasma constituents and /( p) is their 
distribution function. In case of electromagnetic plasma the coupling constant g should be replaced by e but the 
results are then in the so-called Lorentz-Heaviside units not in the Gauss units typically used in the plasma physics. 
For the quark-gluon plasma governed by QCD with the SU(iV c ) gauge group and one quark flavor, /( p) = n(p) + 
n(p) + 2N c n g (p), where n(p), n(p), n g { p) are the distribution functions of quarks, antiquarks and gluons of a single 
color component. The chromodielectric tensor does not carry any color indices, as the state corresponding to the 
momentum distribution /(p) is assumed to be colorless. The z0 + prescription makes the Fourier transformed dielectric 
tensor e lJ (f, r) vanish for t < 0, which is required by causality. In kinetic theory, the infinitesimal quantity i0 + can 
be treated as a remnant of inter-particle collisions. Performing the integration by parts, the dielectric tensor Q can 
be rewritten in the form 


e y (w,k) = 5 ij 


g 2 

f d 3 p 

/( P) 

fcV + Tdfc- 7 (k 2 — 

X l 3 | | v J 

2w 2 J 

(2t r) 3 

E P 

ai — v • k + i0+ ' (ui — v • k + *0+) 2 


( 5 ) 


which will be more convenient for our purposes than the expression Q. 

The dielectric tensor given by Eq. Q or © is fully determined by the momentum distribution of plasma con¬ 
stituents. The distribution function of the two-stream system is chosen to be 


/(P) = ( 27 r) 3 p[< 5 ( 3 ) (p- q) +<i (3) (p + q) 


( 6 ) 
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where p is the effective density of plasma constituents in a single stream. The distribution function ([d]) can be treated 
as an idealization of the two-peak distribution where the particles have momenta close to q or — q but it is not required 
that the momenta are exactly q or q. 

The distribution function © substituted into Eq. © provides the dielectric tensor which in turn gives the matrix 
© of the form 


E y (w, k) = (w 2 - k 2 - p 2 )S ij + Wk? - 


UJ‘ 


M 2 (k-u) 
- (k-u) 


+ u l k 3 ) — 


p 2 {oj 2 + (k • u) 2 )(k 2 - oj 2 ) 
(w 2 - (k ■ u ) 2 ) 2 1 


( 7 ) 


where p 2 = g 2 p/Eq is a parameter analogous to the Debye mass squared, and u = q /Eq is the stream velocity. It is 
smaller than the speed of light, if plasma constituents have non-vanishing mass. If we consider a system of massless 
constituents, Eq = |q| and u 2 = 1. However, when the distribution © is treated as an approximation of the two-peak 
structure and the particles have non-zero momenta perpendicular to the stream velocity, Eq < |q| and u 2 < 1 . We 
also note that p is the only dimensional parameter which enters the problem. Therefore, all dimensional quantities 
can be expressed in the units set by the appropriate power of the mass p. 

Our objective is to find a complete set of solutions of the dispersion equation (|3l with the matrix E given by Eq. 0 
It is not difficult to compute the determinant and solve the dispersion equation 15 ]) for specifically chosen orientations 
of the wave vector. In this way we will consider two special cases in Sec. IV In general, however, it appears much 
easier to invert the matrix S and look for poles of E _1 than to find zeros of the determinant of E. Therefore, we 
invert the matrix E in the subsequent section. 


III. POLES OF S ” 1 

To invert the matrix E, one expresses the matrix together with the inverse matrix E -1 in a given basis. Then, the 
inverse matrix is found directly from the equation EE^ 1 = 1. To choose the appropriate basis, one observes that 
the matrix E given by Eq. © is symmetric (E 1 - 3 = E- 3 *) and that it depends on two vectors k and u. Therefore, the 
matrix can be completely decomposed in the basis of four projectors built out of the vectors k and u and the unit 
matrix. Following Romatschke and Strickland [5], we introduce the vector u t which is a component of u transverse 
to k that is 


k l k 3 , 


4 = (■*« - -jjr) 


and we define four projectors 


PM Uihj 

A i3 = 6 ij -^, rr 3 = * . 

K K 


C ij = 


UirpUrp 


,2 ’ 


D 13 = k l u J T + k°u l T . 


The matrix E is expressed in this projector basis as 

S'-' = a A' 1 + b B ij + c ("■' + dD ij 


with the coefficients a, 6, 
= 

With the matrix E given 
a(co, 

b(uj, 

c(w, 


c, d provided by the equations 

k 2 b, = u|(a + c), = u^k 2 d, TrE = 2 a + b + c. 

by Eq. 0: the coefficients a, 6, c, d equal 
k) = uj 2 — p 2 — k 2 , 

^ _ , ,2_ 2 2/i 2 (k ■ u) 2 p 2 (id 2 + (k • u) 2 ) (k 2 — a; 2 ) (k ■ u) 2 

M w 2 — (k ■ u) 2 (w 2 — (k • u) 2 ) 2 k2 

p 2 {u> 2 + (k ■ u) 2 )(k 2 — ui 2 ) ( (k- u) 2 \ 

(w 2 - (k ■ u) 2 ) 2 ^ k 2 ) ’ 

p 2 (k • u) A« 2 (w 2 + (k • u) 2 ) (k 2 - w 2 )(k • u) 
w 2 -(k-u) 2 k 2 (w 2 - (k-u) 2 ) 2 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


d(ui, k) 


(15) 
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When compared to the matrix decomposition applied in the studies [S] and [5], there is one important difference: 
the vector u, which is analogous to the vector n from [61 and |8], is not of unit length but u 2 < 1. Therefore, the 
expression u 2 shows up in the formula (14). 

Expressing the inverse matrix £ _1 in the same basis © and solving the equation ££ 1 = 1, one finds the inverse 
matrix in the following form 


(y -iyj = 1 A ij , ~ a (° + c) Blj + (~d 2 k 2 u| + be) C ij + adDV 
a a(d 2 k 2 Uy — b(a + c)) 

As seen, the poles of the matrix £ _1 (u;,k) are determined by the equations 

a(w,k) = 0, 

b(u, k)(a(w, k) + c(w, k)) — k 2 u^ci 2 (w, k) = 0, 

which are the two dispersion equations equivalent to the general dispersion equation ([3]). 

The dispersion equation with the coefficient a(w, k) given by the formula ( [l2| ) has the simple solution 

co 2 a (k )=^ 2 + k 2 , 


(16) 


(17) 

(18) 


(19) 


which represents transverse plasmon. A real problem is to solve the dispersion equation (18). 


IV. SPECIAL CASES kiu AND k || u 


Before considering a general solution of the dispersion equation (18), 
k |j u which will help us to analyze the general case. 


we discus here two special cases k J_ u and 


A. kiu 


When 9, which is the angle between k and u, equals 90° that is k • u = 0, the coefficients ( I2][l5 ) simplify to 

a(w,k) = w 2 — i 2 — k 2 , (20) 

b(u, k) = lo 2 — fj, 2 , 


c(w,k) = 


m 2 (w 2 — k 2 )u 2 


d(u, k) = 0. 


The dispersion equation (18) reads 


/ 2 2 \ ( 2 2 j 2 , 

(Cc? — [A — fl — /c 


/i 2 (w 2 — k 2 )u 2 


OJ “ 


= 0 , 


and it has three solutions 


UJ o(k) = 9 2 , u±(k) = ^ (^A 2 + k 2 ± \J( A 2 + fc 2 ) 2 + 4/x 2 u 2 /c 2 

where A = /.<%/! — v?. So, there are three pairs of modes of opposite sign. 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


To clarify a physical character of the solutions (25), we explicitly compute the matrix £ assuming that u = (0,0, u) 
and k = (A;, 0,0). Then, one finds 


£(w,k) = 


oj 2 -fi 2 


0 

- H 2 - k 2 

0 


0 

0 

uj 2 — n 2 — k 2 + 


2 / 2 i 2 \ 

{cj —k 


(26) 


0 
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and its determinant equals 


det[E(w, k)] = (w 2 — /i 2 )(w 2 — /i 2 — fc 2 ) ^w 2 — ^ 2 — k 2 + 


fi 2 (oJ 2 — k 2 )u 2 


(27) 


The general dispersion equation that is detE = 0, gives, as expected, the solutions ( f25j ) and additionally ( [T9| . 
The structure of the matrix ( [26] ) clearly shows that the mode wo(fc) is longitudinal (the electric field is along the wave 
vector) and the remaining modes w a (fc), ui±(k ) are transverse (the electric field is perpendicular to the wave vector). 
The solutions w 2 (fc), Wg(fc) and w 2 (k), which are all positive, correspond to stable real modes while the solution 
w 2 (fc), which is negative, represent two imaginary modes - the Weibel or filamentation unstable and overdamped 
modes. Let us also note that the solutions Wg(fc) and w 2 (fc) cross each other at k = + u 2 . 


B. kllu 


When the decomposition [10]) is used to invert the matrix E, the case k||u needs some care as then Uj- = 0. The 
fact that the vectors k and u are parallel to each other means that the matrix E actually depends on one vector only. 
Indeed, for 0 = 0, when (k • u) 2 = fc 2 w 2 , the matrix equals 


E ij (u, k) = (w 2 - - k 2 )d ij + ^1 - 


2 n 2 u 2 


H 2 u 2 (oj 2 + fc 2 it 2 )( fc 2 — uj 2 ) 


uj 2 — k 2 u 2 k 2 (ui 2 — k 2 u 2 ) 2 

Consequently, one needs only the matrices A and B to fully decompose S i.e. 

E = a A + b B 1 

where the coefficients a and 6, which are found from the equations 

fc* E ij V = k 2 b, TrE = 2 a + b, 

equal 

a (w, k) = u> 2 — /r 2 — fc 2 

2 2 2fj 1 2 k 2 u 2 h 2 u 2 (uj 2 + k 2 u 2 )(k 2 — uj 2 ) 


k l k?. 


(28) 

(29) 


b(u), k) = w 2 — /j. — 


The inverse matrix is 


to 2 — k 2 u 2 


= - A 


(w 2 — k 2 u 2 ) 2 


B. 


(30) 

(31) 

(32) 


We again have the dispersion equation a(w,k) = 0, which gives the solution (19), and instead of the dispersion 
equation (18) we have 6(w,k) = 0 which gives longitudinal modes, as B projects on the direction parallel to k. 

With the coefficient b given by the formula the dispersion equation is 

w 2 (w 4 - (2 k 2 u 2 + m 2 (1 - u 2 ))uj 2 + k 4 u 4 - n 2 k 2 u 2 ( 1 - m 2 )) = 0, (33) 

which, except the trivial solution w 2 = 0, has two solutions 

to±{k) = k 2 u 2 + ± ^\/8 k 2 u 2 + A 2 , (34) 

where, as previously, A = fiy /1 — u 2 . As seen, w 2 (fc) is always positive but w 2 (k) is negative for k 2 < X 2 /u 2 . Then, 
we have the instability which is well known in the plasma physics as the two-stream electrostatic instability. When 
u 2 —» 1, the mode w_(fc) becomes stable and w 2 (k) = w 2 (k) = fc 2 . 

The solutions (191 and (34) can be also easily found directly from the matrix E. Choosing u = (0,0, it) and 
k = (0,0, fc) the matrix equals 


E = 


w 2 - /i 2 — fc 2 0 

0 uj 2 — fjj 2 — fc 2 


0 


0 


2 2 2 fi 2 k 2 u 2 

U ~ H - bJ 2_ k 2 u 2 


2/ 2 i j 2 2 w j 2 2 \ 2 

fl (cj -\-k U )[k —UJ )U 

(uj 2 — k 2 u 2 ) 2 


( 35 ) 


The equation detE = 0 has, as expected, the solutions (19) and (34) and the former one is doubled. 
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FIG. 1: (Color online) The discriminant A as function of k and 8 for u 2 = 1/2 (left panel) and u 2 = 3/2 (right panel). 


V. GENERAL CASE 


We discuss here in full generality the dispersion equation (18 1 with the coefficients a, b , c, d given by Eqs. ( 12p5 |. 
A crucial finding is that the ratio ui 2 / (w 2 — (k ■ u) 2 ) factors out in left-hand-side of the equation. Consequently, we 
get the cubic dispersion equation 


a^x 6 + a 2 X 2 + a\X + ao = 0, 


(36) 


where x = w 2 and the coefficients ao, ai, a 2 , a 3 are the real numbers equal to 

a 0 = -(/i 2 u 2 + (k ■ u) 2 ) (/t 2 (k ■ u) 2 T k 2 ((k ■ u) 2 - /r 2 )), (37) 

ai = (k ■ u) 2 ((k • u) 2 + 2k 2 + /r 2 (l + u 2 )) + /i 2 (l — u 2 )(k 2 + /i 2 ), (38) 

a 2 = —k 2 — 2(k ■ u) 2 + /x 2 (—2 + u 2 ), (39) 

03 = 1. (40) 


As well known, see e.g. [ 9 ], all three roots of a cubic equation can be found algebraically. Since the coefficients 
a 0 , ai, a 2 , 03 are real, the character of the roots depends on a value of the discriminant 


A Q 9 0 Q 0 9 

A = 18aodia 2 a3 — Aa^ao + a x a 2 — — 27 a 0 a 3 . 


(41) 


One distinguishes three cases: 


• if A > 0, the roots are real and distinct; 

• if A = 0, the roots are real and at least two of them coincide; 

• if A < 0, one root is real and the remaining two are complex. 


One shows that the discriminant (411 computed with the coefficients (37p0) is nonnegative for any k, if 0 < u 2 < 1 


but there is a domain of k where A is negative for u 2 > 1. This is demonstrated in Fig. [l] where the discriminant is 
plotted as a function of k and 6 for u 2 = 1/2 (left panel) and u 2 = 3/2 (right panel). Since the stream velocity is 
limited by the speed of light, we always have three real solutions of the dispersion equation ( [36] ). 

The real solutions of the cubic equation can be written down in the Viete’s trigonometric form [Si 


w„( k) = 2\l — cos 


1 /3 VSq\ 27ra 

3 arCC ° S ( 2p§ / ~ ”3“ 


d2 

3a 3 ’ 


(42) 
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FIG. 2: (Color online) Dispersion curves w 2 (k) vs. k 2 at u 2 = 3/4 for four values of the angle 9 equal 0° (a), 30° (b), 60° (c), 
and 90° (d). The red (solid) lines are for w 2 (k), the green (dashed) for tj 2 (k), the blue (dotted) for tu 2 (k), and the orange 
(dashed-dotted) for cu 2 (k). The light cone is represented as a light gray line. 


where n = 1,2,3 and 


P = 


3<Z3ai — aj 
3a§ 


q = 


2a% — 903020.1 + 27a§ao 


274 


(43) 


These formulas assume that p < 0 and that the argument of the arccosine belongs to [—1,1]. These conditions are 
guaranteed as long as A = —a|(4p 3 + 27q 2 ) > 0 which is the case under consideration. 

The complete set of dispersion curves of plasmons predicted by the formulas (19) and (42) is shown in Fig. [ 2 ] for 
u 2 = | and four different orientations of the wave vector k. As seen, the solution w|(k) corresponds to the unstable 
and overdamped modes. We also observe in the figure that the green (dashed) line representing w 2 (k) approaches 
the blue (dotted) line which refers to w|(k). At 9 = 90° the lines hit each other but they do not cross. This is the 
phenomenon of mode coupling which is nicely explained in §64 of the textbook [3] . 

How the general solutions (42) are related to those found in Sec. IV for kf u and k jj u? In the former case, we 
have the relations: 


l (k) = 


i(k) = 


wg(fc) 

for 

k < 

w+O) 

for 

k > 

w+O) 

for 

k < 

Wo(fc) 

for 

k > 


Vl+u 2 ’ 

lIU 

— 777 +^' 


Vl+U 


and wf(k) = w 2 (k). The crossing of the solutions Wnffc) and w. 


3 UJnfl 

limit (k ■ u) —► 0. The general solutions shown in Fig. [2^ do not cross each other. 


— Vi+u 2 ’ 

( k ) derived in Sec. 


IV A 


(44) 


(45) 


actually results from the 


The special solutions obtained in Sec. IV B for 
w|(k) = w+(fc) and w|(k) = w 2 (fc). 


u are related to the general solutions as: w 2 (k) = w 2 (k), 



































1.5 >2 2.0 



FIG. 3: (Color online) Dispersion curves ca 2 (k) vs. k 2 at u 2 = 1 for four values of the angle 6 equal 0° (a), 30° (b), 60° (c), 
and 90° (d). The red (solid) lines are for ca 2 (k), the green (dashed) for w 2 (k), the blue (dotted) for w|( k), and the orange 
(dashed-dotted) for caf(k). The light cone is represented as a light gray line. 


VI. SPECIAL CASES u 2 = 1 


We consider here the case when the stream velocity u equals the speed of light. This case is interesting both 
from physical and mathematical points of view. The dielectric tensor of the two-stream system with u 2 = 1 exactly 
coincides (under the replacement u — > n and 2 fj? —> to 2 ) with that of the plasma with an extremely prolate (infinitely 
elongated in one direction) momentum distribution discussed in the study [ 6 j. Therefore, the plasmon spectra are 
obviously the same. Nevertheless, when the spectrum is found as a limit u 2 —>■ 1, the mode crossing observed in the 
extremely prolate plasma gets a different physical meaning. As explained below, instead of mode crossing we rather 
have the extreme mode coupling mentioned in the previous section. 

The dispersion equation (36) u 2 = 1 is also interesting mathematically. The form of the solutions (|42|) with the 


trigonometric and inverse trigonometric functions is required, if we deal with the so-called casus irreducibilis, when 
three real and distinct roots of a cubic equation cannot be expressed in terms of real radicals. However, a cubic 
equation, which has three real and distinct roots, can be sometimes reduced to a quadratic equation by means of 
the rational root test. Then, all three real roots of the cubic equation are expressed by real radicals and the Viete’s 
trigonometric form is an unnecessary complication. 

Since the coefficient 03 in Eq. (36) equals unity, see Eq. (40), the rational root test suggests to look for a root of 
the equation among the factors of ao given by the formula (37). When u 2 = 1, there is a factor fj, 2 + (k ■ u) 2 which is 
indeed the root of the equation. Consequently, the cubic equation (36) is reduced to a quadratic one, which is easily 
solved, and the three solutions read 


ca 2 (k) = M 2 + (k.u) 2 , (46) 

(k) = - (^k 2 + (k • u) 2 ± ^k 4 + (k • u) 4 + A(i 2 k 2 — 4/x 2 (k • u) 2 — 2fc 2 (k • u) 2 ^j. 


( 47 ) 
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The solutions a>Q(k) and w+(k) are positive for any k and consequently they represent real modes. The modes wo(k) 
and w + (k) cross each other at k = g ■ The solution w 2 (k) is negative for k < n\ tan#| and positive otherwise. It 
represents the Weibel unstable mode and its overdamped partner for sufficiently small wave vectors. When klnor 
9 — 90°, the unstable mode exists for all values of k. When kj|n or 9 = 0° the configuration is cylindrically symmetric 
and there is no instability. This is the situation we have already encountered in Sec EH- the unstable longitudinal 
mode disappears asu-> 1. 


Comparing the solutions (46) and (47) to the general ones given by the formula (421, one realizes that 


J i(k) = 


col(k) = 


<*>o(k) 

for 

k < 


y/2 sin 6 

w+(k) 

for 

k > 


y/2 sin 6 

w+(k) 

for 

k < 


y/2 sin 6 

^o(k) 

for 

k > 


y/2 sin 6 


(48) 


(49) 


and w 2 (k) = w 2 (k). Therefore, one sees that the crossing of the solutions wo(k) and w+(k) is actually an artifact of 
the limit u 2 — i 1. The physical solutions are the combinations (481 and (49). 


The complete spectrum of plasmons, which includes w 2 (k), w|(k), w|(k) and w 2 (k), is shown in Fig. [ 3 ] for four 
different orientations of the wave vector k. The qualitative difference, when compared to the case u 2 < 1, occurs 
when 9 —► 0. Then, the solutions w|(k) and w|(k) merge into one double solution w|(k) = w|(k) = k 2 . 


VII. SUMMARY AND FINAL REMARKS 


We have performed a systematic analysis of plasmons in the two-stream plasma system. The complete spectrum, 
which consists of four pairs of modes of opposite sign, has been found in a closed analytic form. At any orientation 

of the wave vector and u 2 < 1 , three pairs of modes are pure real and one pair is pure imaginary for sufficiently 

small wave vectors. When the wave vector is bigger than the critical one, all four pairs of modes are real. Among the 
imaginary modes, one is unstable and one is overdamped. An interesting feature of the unstable mode is that when 
the orientation of the wave vector changes from parallel to the stream velocity to perpendicular, the mode smoothly 
changes from longitudinal to transverse. In the first case, we have the electrostatic two-stream instability and in 
the latter one the filamentation or Weibel instability. When u 2 —► 1, the electrostatic unstable mode disappears. 
Therefore, when u 2 = 1 and k || u, the whole spectrum of modes is real. 

In the limiting cases k _L u and u 2 = 1, one seems to observe mode crossing. However, a confrontation of the 

special case solutions with the general ones shows that there is no mode crossing but rather extreme mode coupling 
when the dispersion curves hit each other but they do not cross. 

A mathematical structure of the plasmon spectrum is also interesting. In the general case we deal with the casus 
irreducibilis when the three solutions of the cubic dispersion equation cannot be expressed in terms of real radicals. 
Therefore, the solutions are given by the trigonometric and inverse trigonometric functions. However, in the extreme 
cases k _L u, k || u and u 2 = 1, the cubic equation is reduced to the quadratic one and the solutions are expressed in 
terms of real quadratic radicals. 

The two-stream configuration repeatedly occurs in the physics of electromagnetic plasma and thus our results are 
hopefully of some practical interest. In case of the quark-gluon plasma produced in relativistic heavy-ion collisions, 
the two-stream system is rather irrelevant. However, such a system can be treated as a simple but nontrivial model 
of an unstable plasma. We have actually used the model to study the energy-loss of a high-energy parton m and it 
can be also applied to other problems. 
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